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In n — p bilayer systems an exotic phase-coherent state emerges due to Coulomb pairing of n- layer 
electrons with p-layer holes. Unlike Josephson junctions, the order parameter phase may be locked 
by matrix elements of interlayer tunneling in n — p bilayers. Here we show how the phase locking 
phenomenon specifies the response of the electron-hole condensate to interlayer voltages. In the 
absence of an applied magnetic field, the phase is steady-state (locked) at low interlayer voltages, 
V < V c , however the phase increases monotonically with time (is unlocked) at V > V c . The change 
in the system dynamics at V = V c gives rise to a peak in the differential tunneling conductance. 
The peak width V c is proportional to the absolute value of the tunneling matrix element IT12I, but 
its height does not depend on |Ti2 1; thus the peak is sharp for small |Jia|. A sufficiently strong 
in-plane magnetic field reduces considerably the peak height. The present results are in qualitative 
agreement with the zero bias peak behavior that has recently been observed in bilayer quantum Hall 
ferromagnets with spontaneous interlayer phase coherence. 

The idea that in bilayer n—p structures consisting of an electron- conductivity layer (n-layer) and a hole-conductivity 
layer (p-layer) the Coulomb attraction of electrons and holes may lead a formation of electron-hole pairs with spatially 
separated components was put forward rather long ago [1,2]. As a result of Bose-Einstein condensation of these pairs, 
there arises a peculiar superfluid (phase-coherent) state, in which a dissipationless motion of pairs gives rise to equal- 
in-magnitude and oppositely directed supercurrents in n- and p-layers. At present, two variants of the systems have 
been realized experimentally where an excitonic condensate with spatially separated components is formed. In both 
cases, these are closely lying GaAs/AlGaAs double quantum wells, where either interwell excitons are excited by a 
laser pulse [3, 4], or two-dimensional electron layers are formed due to doping. In the last case, the electron layers must 
be placed in a strong magnetic field, normal to the layers, such that the total filling factor should be vt = V\ + v-j = 1 
[5]. Since all these systems have one and the same exciton mechanism for the interlayer phase coherence [6], the 
physical properties of these systems in the coherent state must be qualitatively similar. 

The present paper has mainly been stimulated by recent impressing experiments of Spielman et al. [7, 8] who 
have found that if a bilayer electron system transits into a phase-coherent state (in which the quantum Hall effect is 
observed at Vt = 1), then this transition is accompanied by a sharp rise in the differential tunneling conductance Gt 
at low interlayer voltages V. As the temperature is lowered, this peak of tunneling conductance remains of a finite 
height and width in contrast to the tunneling conductance peak of the Josephson junction. In the parallel to the 
layers magnetic field H the peak height decreases the more drastically, the higher is the field, and at H > 0.6 T the 
peak becomes practically indistinguishable. 

A number of papers have been devoted to theoretical interpretation of the experimental results obtained by Spielman 
et al. For example, Fogler and Wilczek [9] have treated the tunneling conductance peak as a consequence of the 
Josephson effect in a long inhomogeneous junction. In refs. [10, 11], the interpretation of the peak has been based on 
the notion of the finite time of phase coherence. Joglekar and MacDonald [12] have performed both phenomenological 
and microscopic calculations of tunneling conductance Gt value at V = 0. In ref. [13], Gt(V, H) was calculated using 
the phenomenological equation similar to the Landau-Lifshitz equation for the magnetic moment. Such a diversity of 
theoretical approaches in the interpretation of experiment [7, 8] gives impetus to a consistent microscopic consideration 
of the dynamics of phase-coherent bilayer systems, this being just the subject of the present paper. Though we consider 
the n — p system in the absence of a perpendicular magnetic field, the exciton nature of the collective state in all 
the above-mentioned systems encourages us to believe that the present results provide a qualitative description of 
experiments by Spielman et al. [7, 8]. 

An important but still not completely resolved problem for the systems with electron-hole pairing is the problem of 
phase locking by interband transitions [14], which coincide with interlayer tunneling transitions in the systems under 
consideration . The tunneling transitions lift the degeneracy in the phase of the order parameter, thereby locking the 
phase and making it equal to the phase of tunneling matrix elements. The last statement is valid in the absence of 
a parallel to the layers magnetic field. Kulik and one of the present authors [15] have shown that in the magnetic 
field parallel to the layers the phase is locked only in the fields H < H c \. (The critical field H c \ oc |Xi2 1 1//2 , where 
2ia(= \Ti2\e lx ) is the matrix element of interlayer tunneling). At H > H c \ the phase locking is lifted and the phase 
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monotonically changes in the direction normal to the field, this giving rise to spatial oscillations of tunneling current 
(vortex state). The phase locking phenomenon appears to exert an essential effect not only on the thermodynamic 
properties of n — p systems, but also on their kinetics. 

Relying on the microscopic approach, the present paper deals with the response of phase-coherent n — p system to 
the interlayer voltage V. We demonstrate that similarly to the existence of the critical field H c i, in the case under 
consideration there exists the threshold voltage V c (<x |Ti2 1 ) that quantitatively characterizes the degree of phase 
locking in the n — p system. At low voltages, V < V c , the order parameter phase is locked (steady-state), and the 
direct tunneling current is proportional to V. The Ohmic character of a spatially uniform tunneling current at V < V c 
means that in the phase-coherent n — p system there is no dc Josephson effect [16]. (The absence of the dc Josephson 
effect in the two- layer electron system has been established by Joglekar and MacDonald [12]). At voltages V > V c , 
the phase monotonically changes with time, and this results in the tunneling current oscillations with frequency 
u) = e^/V 2 — V 2 (here e is the elementary charge, % = 1) So, at V > V c the n — p system retains the essential feature 
of the ac Josephson effect in superconductors, namely, the tunneling current oscillations at a constant applied voltage. 
At the same time, the dissipative character of the oscillating tunneling current (see below), the nonuniversality of the 
voltage dependence of to, and the presence of threshold voltage V c are specific to phase-coherent bilayer n — p systems. 

Further on, we show that the above-described "liberation" of the order parameter phase at V = V c results in a 
sharp peak of Gt(V), the height of which is independent of | Ti 2 1 , and its width is equal to 2V C , i.e., for small |Xi 2 | 
the peak will be high and sharp. Thus, in our opinion, the nature of the tunneling conductance peak observed in 
the experiments of Spielman et al. is closely connected with the phenomenon of order parameter phase locking by 
tunneling transitions. The experimentally observed suppression of Gt(V) peak with an increasing parallel magnetic 
field [8] also lends support in favor of this interpretation, because, as indicated above, a sufficiently strong in-plane 
magnetic field eliminates the phase locking. 

We are now coming to the analysis of the dynamics of a phase-coherent n — p system in the limit of a high pair 
density, when the average distance between the electron-hole pairs is small as compared to the characteristic pair 
size. The advantage of the high-density limit lies in the possibility of considering the phase-coherent system dynamics 
in the gapless state, when the gap in the excitation spectrum becomes zero under the action of strong depairing, 
and the order parameter A is reduced but remains non-zero [17]. For the n—p bilayer the order parameter is 
proportional to the average (ipi(r, t)^"(r, t)), where ip^(r,t) is the electron creation operator in the layer i. An 
essential simplification consists in the fact that the absence of the gap makes it possible to describe the dynamics of 
the phase-coherent system only in terms of the complex order parameter (A = | A|e ze ) without involving the dynamics 
of the quasiparticle distribution function. 

The dynamic equation for the order parameter of the n — p system was derived by the Green function technique in 
our paper [18] and has the following form: 



- (A — teVA) + {A- B\A\ 2 + D[%- + -(A 1 - A 2 )] 2 }A + ^ = 0. (1) 

or c Qt 

The equation obtained is in a perfect agreement with the general theory of a relaxation of an order parameter near 
a point of a phase transition of the second kind (see, for instance [19]). In accordance with this theory a state of a 
physical system under a phase transition of the second kind can be described by an order parameter, that is nonzero 
below the transition point and equal to zero above this point. An equilibrium value of the order parameter can be 
found from the condition that a variation of a corresponding thermodynamic potential is equal to zero. In the absence 
of the interband hybridization the thermodynamic potential for the condensate of electron-hole pairs with spatially 
separated components can be presented in the form 

F = J \d\[-^ v ~ - c (Ax - A 2 )]A| 2 + A|A| 2 + is|A| 4 |dr (2) 

Expression J2J is similar to the thermodynamic potential for Cooper pairs in the Ginzburg-Landau theory, but here 
the term 2eA is replaced by the term e(Aj — A2). Such a modification is quite natural. Indeed, for the case of electron- 
hole pairs with spatially separated components an electron in the layer 1 "sees" the vector potential Ai, while a hole 
in the layer 2 " sees" the vector potential A 2 . Since the signs of the electron and hole charges are different the vector 
potentials Ai and A 2 are subtracted from each other in eq. In the equilibrium the order parameter A(r) is 

found from the condition 6F/6 A(r) = 0. At small deviation from the equilibrium, when the derivative 5F/S A(r) is 
nonzero but a small one, the order parameter relaxation rate (the derivative dA/dt) is also small. In the mean field 
approximation these two derivatives should be proportional to each other. But it is necessary to take into account 
that due to the gauge invariance of the theory the derivative d/dt can enter into the equation in a combination with 
the term ie(Vi — V 2 ), where V\ and V% are the electrostatic potentials in the layer 1 and 2, correspondingly. As a 
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result in the absence of the interband hybridization one arrives to the equation QJ, where T\% = 0. In the presence 
of the interband hybridization the Hamiltonian of the system contains the terms linear in the order parameter A and 
in the matrix elements T\i and the corresponding conjugate terms (and it means that the thermodynamic potential 
contains the same terms). These terms play the role of a source of electron-hole pairs. They are analogous to the 
terms that appear in the Hamiltonian of a ferromagnet in an external magnetic field. For the case of a magnet it 
results in an appearance of a term linear in the magnetic field in the equation for the order parameter. Since for the 
system considered the matrix element T\2 is analogous to the magnetic field, a term linear in Tyi should appear in 
the equation for the order parameter in the presence of the interband hybridization. We see that eq. contains this 
term, indeed. The microscopic analysis shows that in spite of the phenomenological arguments presented look quite 
general, in reality, eq. Q is valid only in a rather narrow interval of the impurity concentration in similarity with the 
Gor'kov-Eliashberg equation for the superconductors with paramagnetic impurities [20]. 

In the gapless situation under consideration, the coefficients of the dynamic Ginzburg-Landau equation Q have 
the forms A(T) = (2tt 2 /3)t(T c 2 - T 2 ), B = Amr/ZM, D = Pqt/M 2 [18]. Here r is the electron elastic scattering 
time ( for simplicity, it is considered equal to the hole elastic scattering time), T is the temperature {Ub = 1), T c is 
the critical temperature, M = mi + 7712 is the pair mass, m = 77117712/M is the reduced mass of pair, po is the Fermi 
momentum of electrons and holes, ( is the dimensionless constant of the Coulomb interaction [2]. It should be noted 
that eq. (JTJ is derivated by expanding of the anomalous Green function as a power series in (A/T c ) [18]. Since a term 
linear in the matrix element J12 appears in the expression for the order parameter, it is necessary that |Xj_2 1 *C T c for 
the validity of eq. . 

At low fields and currents, the modulus of the order parameter varies only slightly in space and time. Assuming 
|A| to be constant equal to A , the imaginary part of eq. Q can be written as follows 



^-D^-(^-^[Hn])-eV + eV c smcj> = 0. (3) 
Or or <3> 

Here, the gradient-invariant phase = # — x — (2Trd/&o)A z is introduced, d is the interlayer distance, $0 = hc/e is 
the magnetic flux. The unit vector n = (0, 0, 1) is normal to the layers and is directed from layer 1 (electron layer) to 
layer 2 (hole layer). The threshold voltage V c = |Ti 2 |/(eCrAo). 

It is readily seen that in the uniform case eq. (J3J) for the phase is different from the equation <f> = eV that appears 
in a number of papers and is treated as the Josephson relation for phase-coherent bilayer systems. The occurrence of 
the term proportional to \Ti 2 \ in the dynamic equation for the phase radically changes the solutions of this equation. 
Thus, in the absence of external fields the stable steady-state and uniform solution of eq. © is = 0, i.e., 6 = x, 
and this means that the interlayer tunneling transitions hold the order parameter phase locked. Below, we consider 
in detail how the phase locking phenomenon influences the dynamic properties of n — p systems. 

We start from the analysis of dynamics of the n — p system in the phase-coherent state for the spatially uniform case 
in the absence of the magnetic field. Let the n — p tunneling junction is incorporated into the electrical circuit having 
the resistance R and the voltage source £ . The resulting voltage V across the n — p tunneling junction determines 
the difference of electrochemical potentials of the layers and thereby dictates the carrier density in n— and p— layers. 
If Sn is the deviation of the electron density from the equilibrium one, then the equality eV — —5n/N*(0) is valid, 

where the renormalized density of states on the Fermi surface is iV*(0) = 7V(0)(1 + ^^) _1 (C is the capacity of the 
bilayer system related to the unit area; N(0) — m/ix). 

In the approximation linear in T\2 the density of tunneling current from layer 1 to layer 2 is equal to J = J c sin 0, 
where J c = 4eA r (0)|Ti2|Ao/C [18]. The charge balance equation can be written as 



eSSn-\ I c sin = 0, (4) 

R 

where S is the area of the n — p junction, and I c = SJ C . Though below we assume £ = const, it should be noted that 
eq. also holds for the time-dependent voltage source. 

Making use of the relationship between Sn and V, from eqs. (© and Q one can derive the second-order equation 
for the phase <fi. In terms of dimensionless variables, this equation takes on the following form, well known in the 
theory of Josephson junctions 

4> + — 7=(1 + e cos 0)0 + sin0 = p. (5) 
VP 

Here, the following dimensionless parameters are introduced: (3 — e£ c to, e = eV c to, p — £/£ c , where £ c = V c + I C R, 
and to = e 2 Nt,(0)RS . The time is measured in the units, where coq = (efc/io) 1 / 2 • 
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Despite the coincidence of eq.lJSJ with the dynamic equation for the phase difference across the Josephson junction, 
the different meaning of the parameter p entering into these equations leads (as it will be seen from what follows) to 
a substantially different behavior of n — p systems and Josephson junctions. 

A detailed analysis of dynamic states of the system described by eq. © was performed by Belykh et al. [21] Not 
going into details of that analysis, we shall mention its main results. For each value of the parameter e one can find 
the corresponding number f3i. At j3 > /3i (large resistances R), the range of p values is split into three adjacent 
intervals: < p < p c , p c < p < 1 , p > 1 (Pc(P, e) is the bifurcation value of the parameter p [21]). In the first interval, 
there is only one stable solution 4> = arcsin p; in the third interval the only stable state is the limiting cycle embracing 
the phase cylinder. In the intermediate (second) interval the both solutions, (f> — arcsin p and the limiting cycle, are 
stable. This non-uniqueness of the solution of eq. (JSJ results in the hysteresis of current-voltage characteristic (CVC) 
at /3 > /?x- For j3 < fa (low resistances R) the stable solutions will be <fi — arcsin p at < p < 1 and the limiting 
cycle at p > 1, while the interval p with two stable states drops out. Correspondingly, at (3 < fix the CVCs have no 
hysteresis. 

Further on, we find the CVC and the differential tunneling conductance of the n — p system in a simple, but 
physically rather demonstrative, case R = 0. In this limit, no distinction may be made between V and £, V c and £ c , 
and the dynamics of the system may be analyzed relying on eq. (J3J) (without spatial derivatives). Since in the case 
considered we have f3 < 0i, then the hysteresis of the CVC is absent. 

If the system is spatially uniform and the voltage V does not depend on the time the equation @ can be integrated. 
One can see that for V < V c the equation |J5} has the time independent solution <fi — arcsin V/V c . In such a case 
the tunnel current, that does not depend on time as well, is equal to I c smcf>o = I C V/V C = V/R c . This current is 
proportional to the voltage V applied and it is a usual dissipative current. 

The corresponding tunneling conductance is given by 



Gt = dV = R ^ = 4e2N ^ rA o S - 

Note that at V < V c the tunneling conductance is constant and is independent of the tunneling matrix element |Ti2| 
value. This independence of tunneling conductance from | T\ 2 1 , and also its proportionality to Ag(T) are in agreement 
with the result of Joglekar and MacDonald [12] for Gt at V = 0. 
In case of V > V c the integration yields the tunnel current equals to 



tan ^ 

I(t) = 21, 



1 + tan 



2 0(t) 



(7) 



where 



tan 



v c /v + ^i~(v c /vy 



tan 



■-(v 2 -v c 2 )Ht-t ) 



(8) 



One can see that the interlayer current oscillates with the frequency u> = e(V 2 — V 2 ) 1 ' 2 and this current is not a 
sinusoidal one. Due to non-sinusoidal character of the oscillations the average value of the tunnel current is nonzero. 
The average current is the function on the voltage V. 



I=(I C /V C )(V-^V 2 -V 2 ). (9) 

The behavior of the system considered is similar to the behavior of a Josephson junction between two supercon- 
ductors in a circuit, where the junction is connected in series with a resistor and a voltage generator. But in the case 
considered the essential difference is that the resistor (with R c = V c /I c ) is embedded in the junction and it cannot be 
deleted from the circuit. There is not any transverse superconductivity in the systems considered. 

Since according to eq. © the tunneling current decreases with an increasing voltage, the differential tunneling 
conductance at V > V c is negative: 



G T (V) = -(I C /V C )[V(V 2 - v 2 )- 1 ' 2 - 1]. (10) 

The conductance Gt(V) has its maximum (constant) value at |V| < V c and the discontinuity points at V = ±V C . 
At |V| > V c , as |V| increases, the tunneling conductance monotonically tends to zero, remaining negative. If we take 
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into account the fluctuation smoothing of the CVC, then the dependence of Gt on V will look like a smooth curve 
with the maximum at V = (approximately 2V C in width) and two minima at V w ±V C . It is just this behavior of 
the Gt(V) curve that was observed in experiment [8] in the absence of the magnetic field parallel to the layers. 

It should be noted that both at V > V c and V < V c the spatially uniform tunneling current is dissipative. The 
reason for the dissipation lies in the fact that the uniform interlayer current causes the order parameter phase to 
deviate from its equilibrium value and a continuous input of energy is required to maintain this nonequilibrium state. 

Let now the bilayer n — p structure be placed in a magnetic field H parallel to the layers and directed along the 
x-axis. If H > H c i = {2^ a /ir 2 d)(J c M/en s ) 1 l 2 (the two-dimensional density of the pairs n s = 4pg7V(0)(rA ) 2 /M), 
then the magnetic field between the layers has a nonuniform (vortex) component. We shall show that the CVCs of 
the n — p system in the magnetic field strongly differ from the CVCs in the zero field and are substantially different 
at both low and high resistances R. In the limiting case R = (and H ^> H c i), the solution of eq. © can be derived 
using the perturbation theory. Putting <p = fa + fa, where fa = ky + ut ( k = 2irdH/<&Q, lu = eV) and taking 
into account the correction term fa (proportional to a small T\i value) as a perturbation, we obtain the following 
expression for the average tunneling current density: 



eV c uj 
J = Jc ~(Dk 2 Y+uj 2 - 

So, for R = the CVC has a wide diffusion maximum at u) = Dk 2 . 

At high R values, the charge transport from one layer to the other over the electrical circuit is insignificant, and 
the electron density dynamics in layer 1 is determined by the continuity equation 



e5n = div2j + J c sin<^, (12) 

where div2j denotes the two-dimensional divergence of the intralayer current j = — ^ff (fjr — j^[Hn]). In the assump- 
tion that eV <C tAq, the above-described perturbation-theory procedure yields the following equation for fa: 



v r,d 2 fa 2 d 2 fa J c 

* D ^ u o W - -^uo) sm<t>0 ' (13) 

where u = (n s /MN r (0))^ 2 . Unlike the R = case, the left-hand part of eq. I|13(l has a wave character rather than 
a diffusion character. Correspondingly, the expression for the average tunneling current density 



J = J r 



2X 2 (u 2 /u 2 



k 2 ) 2 + (tuak 2 ) 2 



(14) 



has the resonance at w = uok, the width of which is determined by the attenuation a = D/v,q. This resonance results 
from the coincidence between the plasmon velocity uq in the bilayer structure and the velocity of the magnetic-field 
vortices. The parameter \j equals (eris/MJc) 1 ^ 2 . 

From relations (|llfl and (|14fl it follows that at H 3> H c i the Gt(0) value is proportional to a small \Ti2\ 2 value, 
i.e., the differential tunneling conductance peak (occurring at H = 0) is strongly suppressed. The reason for this 
suppression lies in the fact that at H > H c i the phase <f> monotonically varies with the coordinate, and in this case 
eq. has no stationary solution at finite voltage, i.e., no phase locking arises. 

Thus, the present work has demonstrated in the frame of the consistent microscopic approach that in phase-coherent 
bilayer n — p systems the known phenomenon of order parameter phase locking by matrix elements of tunneling T±2 
leads to a sharp peak in differential tunneling conductance Gt(V) at V = 0. The peak height is independent of |Ti2 1 
and its width is proportional |2i2|, i.e., at weak tunneling the peak is high and sharp. These results are in qualitative 
agreement with the peculiarities of Gt(V) observed in electron bilayer systems in the regime of an integral quantum 
Hall effect at the total filling factor vt = 1. We stress once again that though the theory developed here describes the 
ri — p system without a transverse magnetic field, the present results are in qualitative agreement with the data from 
experiments on electron bilayer systems in a strong transverse magnetic field. This agreement does not seem to be 
accidental. The reason is that most likely the strong magnetic field does not affect the structure of the equation that 
defines the dynamics of the order parameter, but only changes the values of coefficients entering into this equation. 
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